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Abstract. 

Baryon Acoustic Oscillations (BAO) in the radial direction offer a method 
to directly measure the Universe expansion history, and to set limits to space 
curvature when combined to the angular BAO signal. In addition to spectroscopic 
surveys, radial BAO might be measured from accurate enough photometric 
redshifts obtained with narrow-band filters. We explore the requirements for 
a photometric survey using Luminous Red Galaxies (LRG) to competitively 
measure the radial BAO signal and discuss the possible systematic errors of 
this approach. If LRG were a highly homogeneous population, we show that 
the photo-z accuracy would not substantially improve by increasing the number 
of filters beyond ~ 10, except for a small fraction of the sources detected at 
high signal-to-noise, and broad-band filters would suffice to achieve the target 
(Tz = 0.003(1 -|- z) for measuring radial BAO. Using the LRG spectra obtained 
from SDSS, we find that the spectral variability of LRG substantially worsens the 
achievable photometric redshift errors, and that the optimal system consists of 
~ 30 filters of width AA/A ~ 0.02. A S/N > 20 is generally necessary at the 
filters on the red side of the Ha break to reach the target photometric accuracy. 
We estimate that a 5- year survey in a dedicated telescope with etendue in excess 
of 60 deg^ would be necessary to obtain a high enough density of galaxies to 
measure radial BAO with sufficiently low shot noise up to 2; = 0.85. We conclude 
that spectroscopic surveys have a superior performance than photometric ones for 
measuring BAO in the radial direction. 



1. Introduction 

An important observable to constrain the nature of dark energy is the Hubble 
parameter H{z) [T], since it constitutes a more direct probe to the dark energy equation 
of state than the angular diameter distance da{z) or the luminosity distance d^^z), 
which depend on an integral of H(z). Einstein's equations imply that a homogeneous 
and isotropic universe, which is described by the FRW metric, that is composed of 
matter and dark energy with equation of state pQ = wq{z)pq expands according to 
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where the subscripts Q, fc, M and T refer to dark energy, space curvature, matter, and 
total energy density, respectively. In the absence of space curvature, the quantities 
dA{z) and ^^(z) are related to H{z) via dA{z){l + z) = dL{z)/{l + z) = dz' /H{z'). 

The acoustic oscillations in the photon-baryon plasma, observed as acoustic peaks 
in the CMB power spectrum, are also imprinted in the matter distribution at the 
scale of the sound horizon at the radiation drag epoch, when baryons were released 
from the photon pressure. The mass distribution is traced by the distribution of 
galaxies, and the Baryon Acoustic Oscillations (BAO) can be observed as a peak in 
the galaxy correlation function or as a series of harmonic oscillations in the galaxy 
power spectrum. The sound horizon at the radiation drag epoch can be computed very 
accurately from CMB observations (e.g., = 153.2 ±2.0 Mpc from WMAP5, [Ij), so 
it provides a natural standard ruler. In fact, the galaxy power spectrum can be used 
to measure both the angular diameter distance through the clustering perpendicular 
to the line-of-sight, and the expansion rate H{z) through the clustering along the line- 
of-sight. Therefore, BAO measurements test the relation between H{z) and dA{z), 
providing constraints on dark energy and a limit to space curvature (e.g., [3J|31[51|B] ). 
The BAO technique is being considered a powerful probe to the nature of dark energy 
[3 m 13 HO] because of its potential to provide a standard ruler at different redshifts 
and its robustness to systematic effects. 

For an ideal galaxy survey, [TT] have shown that a volume of 1 (Gpc/h)'^ at low 
redshift can constrain Hq to the 7 % level. Forthcoming surveys with larger volume 
are expected to reach the statistical power to constrain H{z) at the % level. There 
are a number of requirements that these surveys need to satisfy for measuring BAO: 
covering a large survey volume, modeling the effects of galaxy bias and non-linearity, 
characterizing the covariance between different modes, evaluating the galaxy selection 
function to sufficient accuracy, reducing photometric calibration errors to low enough 
levels to avoid contamination of the BAO signal, etc. Measuring the BAO scale 
in the radial direction demands in addition that galaxy redshifts are measured to a 
sufficiently high accuracy, ct^, to avoid an excessive smoothing of the BAO peak, which 
has an intrinsic width 5r ~ 10 Mpc, i.e., 5z < SrH{z)/c. As shown by |T2j, a redshift 
accuracy < 0.003(1 + z) is required to avoid substantial loss of accuracy of the 
H(z) measurement from a given survey volume. Spectroscopic surveys usually yield 
a redshift accuracy much higher than this minimum requirement. 

An alternative approach to measure the large number of redshifts required 
for BAO detection are photometric redshifts from imaging surveys. Broad-band 
photometry with ~ 6 filters usually reaches only to ^^/(l + z) ^ 0.03 for the general 
population (with red galaxies having slightly smaller photo-z errors than blue ones), 
insufficient for measuring radial BAO. However, as the Combo-17 survey |13j has 
demonstrated, the galaxy photo-z accuracy can be improved by using a larger number 
of narrower filters. Photometric surveys can cover a large area of the sky faster 
than a spectroscopic survey and reach a higher number density of observed objects. 
We are therefore motivated to investigate the requirements for a photometric survey 
with medium to narrow bands to deliver interesting BAO measurements, and the 
optimization of the number of filters. Previous work has already explored this issue 
f |14l [T2l [151 [T6|). Here we concentrate specifically on the impact of the number of 
bands, signal-to-noise, and non-uniformity of the galaxy sample. We conclude with 




Photo-z optimizaton 



3 



several considerations on the systematic effects that the photometric approach entails. 

In our investigation we use both synthetic stellar population models and Sloan 
Digital Sky Survey (SDSS)-DR6 spectra of Luminous Red Galaxies. We concentrate 
on LRG because they have several properties that make them particularly useful for 
BAO surveys: they are a fairly homogeneous population, the form of their spectra 
makes them particularly suitable for good photo-z determinations, and their high 
luminosity facilitates reaching a high enough signal-to-noise up to high redshifts. For 
these reasons, LRG are the target of choice for z < 2 BAO surveys (at higher redshift 
the Lya forest probably provides the best method for measuring BAO; see [T7]). The 
main conclusion we reach is that a spectroscopic survey is superior to narrow-band 
photometric surveys for measuring the radial BAO signal. We also show that if it 
were possible to find a very homogeneous population of LRG (with spectra closely 
matched by a single spectral template) , then the photometric redshift accuracy would 
not substantially increase with the number of filters used beyond a total number of 
~ 10 at a fixed total exposure time. This is not the case for galaxies measured at high 
signal-to-noise, for which a larger number of filters is optimal, but this high signal- 
to-noise cannot be achieved for a large enough number of objects in a way that is 
competitive with the spectroscopic approach. In reality, however, the variability of 
realistic galaxy spectra worsens the photometric redshift accuracy, making it optimal 
to increase the number of filters to ~ 30 and requiring a higher signal-to-noise per filter 
to reach the desired redshift accuracy. These results are generally in good agreement 
with those of [15] . 

This paper is organized as follows: in § 2 we review the requirements for a survey 
to measure radial BAO. The modeling of the LRG population is described in § 3, and 
in § 4 we describe our fiducial survey model. The results are presented in § 5, where 
we analyze in detail the photo-z accuracy as a function of the number of filters, galaxy 
luminosity and redshift, first for ideal galaxies that match the templates precisely and 
then for real galaxies with SDSS spectra. Discussion and conclusions are presented in 
§ 6 and § 7. Throughout this paper, we use a cosmological model with Hq — 70 km s~^ 
Mpc~^, ^Im = 0.3, r^A = 0.7. Readers who want to quickly see the main conclusions 
or our study may wish to go directly to Figure 14. This shows the number density of 
LRG in several redshift bins with photometric redshift better than 0.003(1 + z) as a 
function of the etendue times the exposure time of a survey. It also shows the number 
density required to reach nP = 1, necessary to make shot noise subdominant in the 
Fourier modes near the line-of-sight useful for measuring BAO. 

2. Spectroscopy vs photometry and target photometric requirements 

Future BAO surveys need to: a) cover large volumes of the universe sampling 
the acoustic scale — 153.2 ± 2 Mpc with a density of galaxies high enough to 
make shot-noise subdominant, and b) not degrade with redshift errors the line- 
of-sight information that yields a measurement of H{z). The latter condition 
implies, quantitatively, that the statistical photo-z errors need to be smaller than 
(Tz — 0.003(1 -t- z) [12]. In Figure [l] we show the present-day (z=0) linear matter 
power spectrum (dotted line) and the corresponding line-of-sight redshift-space power 
spectrum (thin solid line), assuming a scale-independent bias b — 2, a /3 parameter 
f3 = 1/bdlnS/dlna ^ /b ~ 0.24, and assuming an intrinsic galaxy velocity 
dispersion of 420 km/s, corresponding to a pairwise velocity dispersion ap = 600km/s, 
typical of non-linear redshift-space distortions. The line-of-sight power spectrum 
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Figure 1. Galaxy power spectrum in the LCDM model. Dotted line: Present- 
day, real-space, linear matter power spectrum. Solid lines: Redshift space galaxy 
power spectrum on the line-of-sight, for a scale independent bias 6 = 2 and 
redshift error az = 0.003 {thick line) or assuming typical small-scale redshift space 
distortions only {thin line). Long dashed lines: Angle-averaged power spectrum, 
also including redshift errors {thick) or not thin. The shot-noise contribution level 
for a galaxy density of 5 X 10~^h^ Mpc~^ is shown by the short-dashed horizontal 
lines for reference. 



is also shown with Gaussian photo-z errors of CTz ~ 0.003(1 + z) added {thick 
solid lines). The effects of biasing, linear redshift-space distortions [18 , and the 
combined dispersion from photo-z errors and non-linear redshift-space distortions can 
be modelled as: 

Pg(fc, Ai) = 62(1 + ^?l3fP{k)D{ka,ti) , (2) 

where ^ = cosfl, and 9 is the angle with respect to the line-of-sight. In equation [2] the 
"Kaiser factor" is strictly valid only for linear, large-scale redshift-space distortions, 
but it will be sufficient for our main purpose of illustrating the requirements on the 
redshift accuracy and the shot-noise that need to be reached. The function D{kazy) 
describes the small-scale smearing of power. It is often parameterized by a Gaussian, 
D{k(TzfJ,) — exp[— (fccr^/i)^], or by a Lorentzian, D{kazH) = [1 + {ka^n)'^]"^ . If 
X = ka^n << 1, the two descriptions are equivalent since exp(— a;^) ~ (1 + x'^)^^. 
The effect of small-scale nonlinear redshift-space distortions is approximated as a 
contribution to CTz (subdominant in this case compared to redshift errors) . Assuming 
fccr^/x ^ 1 and using the Lorentzian form for D[kazH), the angle-averaged (monopole) 
power spectrum can be calculated analytically with the following result: 

J ^Pg{k,fl)dfi~ (3) 
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This monopole term is shown as the long-dashed line in Figure [T] The Figure 
also shows the shot noise contribution (short-dashed lines) for a galaxy density of 
n = 5 X 10""^, and the level at which hP — 5 and hP = 1 , which is where the 
shot noise contribution is negligible: in the line-of-sight, shot noise starts becoming 
comparable to the P{k) signal at larger scales than for the angle-averaged P{k). 

Let us consider the following two cases: a photometric survey with target 
photo-z errors az ~ 0.003 (thick lines), and a spectroscopic survey (thin 
lines, corresponding to unavoidable non- linear velocities). The sampling variance 
error on P(k,^) scales as apjP = 1/ y/2Trk'^AkAiJ,Vef j /(2'k)^ , where Vefj — 
[fiP(k, ii)/(hP(k, n) + 1)]^ Vsurvey, and the accuracy at which the acoustic scale can 
be measured is directly proportional to ap/P. Thus, if the two surveys have the 
same number density of galaxies, and the spectroscopic survey reaches nP — 5 along 
the line of sight at fc = 0.32/i/Mpc, then the photometric survey reaches nP = 5 at 
k — 0.16ft-/ Mpc. Note that the number of independent modes is roughly proportional 
to fcmoa;, and so the spectroscopic survey could obtain a much better constraint on 
the BAO scale from measuring the power on many more modes, [f] In order to be 
competitive with a spectroscopic survey, a photometric survey would need to achieve 
a much higher galaxy density. If the survey volume were to be the same for the two 
surveys, then at A: = 0.2/i/Mpc a photometric survey would need a galaxy number 
density exp(— fc^cr^) ^ 25 times higher than a spectroscopic survey to achieve the same 
nP along the line of sight direction. The angle-averaged quantity is, of course, less 
sensitive to the smearing along the line of sight, and the same is true for all orientations 
where fj, < 1. As shown by [L2^, at z~l a redshift error of 0.3% degrades the error 
on H(z) by a factor ~ 2, demanding therefore a survey with 4 times the volume of 
a spectroscopic survey to match its performance. As most forthcoming spectroscopic 
surveys will cover more than a quarter of the available sky (i.e., more than a quarter of 
the 30000 square degrees outside the galactic plane), the photometric approach might 
only be advantageous if it could reach higher redshifts. 

The above considerations indicate that a spectroscopic survey is the favored 
option unless a much larger fraction of the sky can be covered with a photometric 
survey and with an extremely high object density; this is equivalent to imposing the 
requirement of cr^ < 0.3% down to fainter magnitudes. Below, we calculate if such a 
survey is possible. For this we concentrate on LRG, which are bright and have very 
homogeneous spectra, and we study the dependence of the results on intrinsic galaxy 
variability, luminosity, etc. 



3. Models for the Population of Luminous Red Galaxies 

We start by calculating the number density of LRG that can be observed at each 
luminosity and redshift. We use the jTHI luminosity function of LRG and adopt 
their model of a Shechter luminosity function slope of a = —0.5 (see Table 6 in 
[19]). To model the spectral distribution of galaxies, we use the first five templates 
presented in [20] [|] These are empirical templates computed using the SPEED [2T] 

J Since the relative error on the power spectrum crp/P does not depend on redshift, the above 
considerations are valid at any redshift, with the caveat that the non-linearities decrease with z but 
the photo-z errors do not. 

§ http:/ /www. ice. csic.es/personal/jimenez/PHOTOZ/ 
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Figure 2. Spectral energy distribution of the five galaxy templates used in this 
paper, labelled from to 4, from 1201 . For displaying purposes, in these figures a 
20 A smoothing has been applied to the lines. 



and CB08 (Chariot, private communication) models with solar metallicity and with 
star formation histories that are representative of the observed distribution in red 
galaxies over the stellar mass range 10® to lO^^M© extracted from [22]. The spectral 
templates correspond to the first five curves of Figure 4 in [22] . The dust reddening 
is performed using the method described in [22] . 

The spectra of the five templates, labeled by the numbers to 4, are shown in 
Figure |2]in the wavelength region of interest for this work. Template is the spectrum 
of the oldest population and was build using the SPEED models, which have a different 
physical treatment for the giant and horizontal branches than the CB08 models. The 
average stellar population age is gradually reduced for the other templates, which 
use the CB08 models. We explain in the next section how we simulate the galaxy 
photometry using these spectral templates. We use the radius-luminosity relation 
measured in SDSS [HdH], 

log(re) = 1.0 - 0.26{Mb - 5 log /i + 21.81) , (4) 

where the effective radius rg is expressed in kpc and Mg is the B-band absolute 
magnitude. We include passive luminosity evolution (with fixed radius) as described 
by a [21] model where star formation starts at z = 4 and decays with a timescale 
r = 0.6 Gyr. This model provides a good match to the color evolution of LRG, 
as reported by [19]. A fixed aperture of 2 arcsec^ will be assumed to simulate the 
galaxy photometry. This represents a reasonable compromise between minimizing the 
amount of galaxy light lost outside the aperture, and the noise contributed by the 
sky. Photometric redshifts are also most accurate when measured on fixed apertures, 
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Figure 3. Fraction of light of a LRG within a 2 arcsec^ aperture as a function of 
luminosity for four different values of redshift, z = 0.3, 0.5, 0.7, 0.9 from bottom 
to top, and two different values of the seeing, 0.8 arcsec (solid lines) and 1.2 arcsec 
{dotted lines). 

rather than variable apertures that may be adjusted to the observed galaxy profile. 
We show in Figure [3] the fraction of galaxy light included within the aperture as a 
function of luminosity for four different values of redshift and two different seeings. 

4. Survey Model 

Throughout this paper we consider a fiducial narrow-band photometric survey as an 
example of the accuracy that can be achieved to measure a large number of LRG 
photometric redshifts, for the purpose of detecting BAO. We focus on narrow-band 
photometry in a wavelength range that is most useful for obtaining LRG redshifts over 
the range z = 0.5 to z = 1.0, using the 4000 A (or Ha) break, a blend of H and Ca lines. 
We choose a set of Nf filters covering the fixed wavelength range from Ai 5300 A to 
A„ = 8300 A, dividing this range into Nf intervals of equal wavelength width. We will 
consider different values Nf to optimize this number for photometric redshifts for a 
given survey and telescope setup. The shape of the filter window functions is assumed 
to be a top-hat of width AA = (A„ — Xi)/{Nf + 1/2)[|[| with the addition of lateral 
wings on each side where the window function varies linearly from zero to the value 
in the central top-hat. The width of each of the wings is set to 1/4 of the width of the 
top-hat. The top-hat parts of the window function are adjacent and non-overlapping, 
whereas the wings cause the filters to have a certain degree of overlap. In a practical 

II The 1/2 is due to the presence of the wings. 
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Figure 4. Overall throughput {solid) and throughput without including 
atmospheric absorption (dotted) used in this work. 

application, it would probably be useful to complement this filter system with wider 
filters around the wavelength range considered here: this would eliminate some of the 
photo-z catastrophic failures (outliers in the photometric redshift error distribution). 
However, the addition of wide filters would not improve the accuracy of the good 
redshifts because they cannot contribute to a better resolution. 

Figure |4] shows the overall throughput as a function of wavelength. We have 
included atmospheric absorption assuming an average of 1.2 atmospheric columns. 
We also include two mirror reflections, filter transmission and CCD efficiency from 
|26l l27l [28] . We use the sky brightness of Patat (2008, private communication; [25]). 
measured for the Paranal Observatory, to simulate the photometric noise in each 
band. This is adequate for a dark site in the absence of any artificial light, and a 
mean air glow intensity over the solar cycle. Figure |5] shows the spectrum of the sky in 
the AB magnitude system, and the average sky brightness in each one of the filters for 
the case Nf = 30 {blue diamonds). The squares and circles indicate the noise levels in 
our fiducial survey and are discussed below. 

Once the shape of the filters, the overall throughput, sky brightness, source flux 
and radius, and seeing have been fixed, the photometric precision reached by a survey 
for the flux in each filter is proportional to [Et/{Nfil,)]^^'^ , where E is the etendue 
(the product of the effective telescope aperture times the field of view) , t is the total 
survey observing time, and is the total solid angle covered by the survey. We shall 
assume for our fiducial survey a product Et = 1.5 x lO^m^deg^hr, corresponding 
for example to a characteristic case of a dedicated telescope with etendue E — 30 
(e.g., a 3 m telescope with a 5 square degrees field of view) that can obtain good 
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Figure 5. Solid line: Paranal sky spectrum from 1251 . Diamonds: Sky surface 
brightness in AB mag/arcsec^ in 30 equally-spaced narrow-band filters. Empty 
circles: 5a sky noise for an aperture of 2 arcsec^ in 30 narrow-band filters using 
a constant exposure time in our fiducial survey. Full circles: 5cr sky noise for an 
aperture of 2 arcsec'^ in 30 narrow-band filters using a linearly growing exposure 
time, with five times longer exposure in the reddest filter than in the bluest one. 
Squares: 5a magitude limits, with the same exposure times as for the full circles. 



photometric data for 1000 effective hours of observation per year, over a period of 5 
years. We assume = 8000 square degrees, close to the maximum area observable 
from a non-equatorial site, since, in general, the best strategy for measuring BAO is 
to cover the widest possible area. We also assume an optimistic value for the seeing of 
0.8", corresponding to the best observatories in the world. Figurejsjcan be used to see 
how the signal-to-noise is degraded for larger seeing: the required value of Et to reach 
a fixed signal-to-noise scales as the inverse square of the fraction of light within the 
aperture. For example, a change of the seeing from 0.8" to 1.0" corresponds to a ^ 
10% reduction in the fraction of light within the aperture, which must be compensated 
by an increase of 20% in Et. 

The number of useful observing hours per year obviously depends on the quality 
of the observing site. The total number of night hours per year is ~ 3500; the number 
of observed hours is reduced owing to cloudy or non- photometric (i.e., affected by 
cirri) nights, and Moon time. Ideally, the presence of increased sky brightness (either 
artificial or lunar) simply means that the effective observed time is reduced as the 
inverse of the sky brightness, because the photometric noise is inversely proportional 
to the square root of the ratio of the sky brightness over the exposure time, when the 
noise is dominated by the sky. For example, an increase by half a magnitude of the 
sky brightness needs to be compensated by an increase of Et by a factor 10°'^ ~ 1.6. 



Photo-z optimizaton 



10 



This implies an effective reduction of the observed time due to the Moon of ~ 25%, 
averaged over the Moon cycle, in the wavelength range of interest. Our assumption 
of 1000 effective hours of observation per year therefore assumes that ~ 40% of the 
night time is clear and with acceptable photometric conditions. The need to co-add 
images obtained under different seeing, sky brightness and transparency conditions 
also needs to be taken into account as an added difficulty. In order to do faint galaxy 
photometry, one often needs to discard the images with the worst seeing and convolve 
the rest to a common maximum seeing, in order to avoid systematic photometric errors 
that depend on galaxy morphology and may tend to introduce artificial correlations 
in the photometric rcdshift errors. The assumptions we make here should therefore 
be regarded as optimistic and corresponding to an excellent observing site. 

For our chosen value of Et = 1.5 x 10^, and seeing of 0.8", Figure [5] shows the 
5 — 17 sky noise in AB magnitudes within a fixed aperture of 2 arcsec^. The open 
circles show the sky noise when the total exposure time is divided equally among 
Nf = 30 filters, and the filled circles assume the exposure time is a linearly increasing 
function of wavelength to reduce the noise in the reddest filters and achieve a more 
uniform LRG redshift accuracy over the range 0.5 < z < 0.9. For the latter case, we 
set the exposure time for the reddest filter at 5 times that for the bluest one. We 
have found these distribution of exposure times to be a reasonable compromise to 
achieve the best redshift accuracies for the largest number of galaxies. We shall use 
these variable exposure times for our fiducial survey throughout the paper, with 5-ct 
sky noise indicated by the filled circles. To compute the total noise in our simulated 
photometry, we add quadratically the sky, source and read-out noise. The read-out 
noise in a photometric mesurement depends on the number of exposures, the pixel size 
and the quality of the CCDs. We model the read-out noise as 7 electrons per pixel 
at each exposure, with a pixel size of 0.4 arcsec. We assume that three exposures 
are obtained for each field and filter (in practice the number of exposures could not 
be reduced below 3 in most filters if the telescope operates in passive drift-scanning 
mode, the best strategy to minimize photometric calibration errors). The magnitude 
of a source detected at a signal-to-noise of 5 under these assumptions is shown as the 
filled squares in Figure [5j The squares approach the filled circles (the 5-a sky noise) 
as the read-out and source noise become small compared to the sky noise. 

The 1-a sky noise is shown in Figure [6] as filled circles. The curves show the AB 
magnitudes of an galaxy in each of the Nj — 30 filters, at four different redshifts. 
Note that the true S/N measured for an galaxy is slightly decreased compared 
to that shown in this Figure because of the source photons contribution to the total 
noise; this is only a ^ 10% effect for an galaxy aX z — 0.5 on the first filter to the 
red of the Ha break. 

To model the accuracy of photometric redshifts, we proceed as follows in the next 
section: first, we simulate the iVj-filter photometry of a galaxy that has the spectrum 
of one of the five templates we use. We add noise for a survey with the characteristics 
described above, and fit the resulting magnitudes with the same template. This will 
tell us the minimum photometric redshift error that would be obtained in the ideal 
case of a galaxy for which we know exactly the spectral shape, and we fit only one 
parameter, the redshift, to the measured photometry. We refer to this as "single 
template case" . We shall examine the dependence of the results on the template used, 
and comment on the effects of galaxy evolution. We shall then consider the case where 
a linear combination of the five templates is used to fit to the photometry simulated 
with one template (5-templates case) and the case where the full set of 5 templates is 
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Figure 6. Simulated photometry of L* galaxies using our template 1 at 
z = 0.5,0.7,0.9,1.1 {solid lines), and the Icr sky noise (circles) and magnitude 
limits (squares). 

used to simulate and fit the photometry (5-to-5 templates case). This will show how 
the redshift error increases when galaxies are allowed to have variable spectra with 
five parameters (the redshift and four relative template amplitudes), but the spectra 
are still assumed to be precisely modelled by the templates used in the fit. In these 
two cases, the effect of changing the number of filters will be explored. 

Finally, we shall use observed spectra of LRG from the SDSS public data to 
generate a photometric "mock" catalog, and fit their photo-z's using our 5 templates. 
The resulting photometric errors should be more realistically expected for real galaxies, 
which are not perfectly matched by a limited number of templates. We note that the 
5 templates we use have been optimized to fit the spectra of a wide class of observed 
galaxies, and therefore they should be an ideal set for obtaining the best possible 
photometric redshifts. 

5. Photometric Redshift Accuracy 

5.1. Optimal number of filters: analytic considerations 

Before presenting our numerical results on the simulated photometry described in 
the previous sections, it will be useful to derive some analytical results on the 
photo-z errors and the optimization of the number of filters under some simplifying 
approximations, using simple examples of galaxy spectra. We consider here the ideal 
case where a galaxy has exactly the same spectrum as the template used to fit the 
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measured photometry, and approximate the filter windows as a top-hat. 

If a fixed total exposure time is to be divided among all the filters, then the signal- 
to-noise in each individual filter, F/ap (where F denotes flux and <jf its statistical 
error), is inversely proportional to the number of filters Nf. the exposure time in 

each filter is tj cx NJ^, and the wavelength width of each filter is A A cx NJ^. 
Hence, the number of photons detected in each filter is Nph oc t/AA cx Nj'^, and 

we increase the number of filters, the resolution increases 
at the expense of the achieved signal-to-noise. 

To model the use of the Ha break in LRG, we first consider a galaxy spectrum 
that has a break at wavelength Aq, where the flux per unit wavelength is F\ = Fq at 
A > Ao and = (1 — B)Fo at A < Aq. The filter that includes the break wavelength 
Aq has a flux 

F = Fo(l - Bx) , (5) 

where x = (Ai — Aq) /AA, and Ai is the wavelength of the right edge of the filter. From 
the measured value of F (and assuming that we know exactly the values of Fq and B 
from the measurements in the other filters), the wavelength Aq can be measured to an 
accuracy 

.. = AA^. (6) 

Since ap/Fc cx Nf, and AA cx NJ^, we conclude that the error to which the 
wavelength Aq is measured is independent of the number of filters, and results in 
a photometric redshift accuracy given by 

cr^ _ ax AA ap 

In other words, unresolved breaks in the spectra of galaxies yield a redshift 
measurement that does not improve as the number of filters is increased. This is 
correct only in the idealized situation where the amplitude of the break is perfectly 
known and galaxy variability may be ignored. 

For a resolved break, we assume that the flux varies linearly from Fq to _F'o(l ^ B) 
over a wavelength range Aq to Aq + SX. The flux measured in a filter centered at 
wavelength A^ and which is fully included within the interval of the break of width 
SX is F — Fo{l — Bx), where x = {Xi — Xo)/6X. Hence, the error on Aq that can be 
deduced from the flux measured in one filter only is 

<7x = SX^ocNf. (8) 

Since the number of filters contained in the wavelength range of the resolved break, 
5X, is proportional to Nf, and the combined error from the measurement in all the 

— 1/2 

filters is reduced as Nj , we find that the set of all photometric measurements will 

1 /2 

yield an error ax cx Nj^' . Therefore, the filters should not be any narrower than the 
intrinsic width of any breaks that substantially contribute to the photometric redshift 

measurement. 

On the other hand, if the spectrum contains features similar to an emission or 
absorption line, then it is advantageous to increase the number of filters. For a line 
feature with an equivalent width Wx which is entirely contained in one filter, the 
measured fiux is F = Fo(l + Wx/^X), independently of the position of the line within 
a top-hat filter. The difference F — Fq = FqWa/AA oc Nf can be measured to an 
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accuracy ap oc Nf, so the position of the feature is measured to increasing accuracy 
as the number of filters is increased, up to the point where the feature is resolved. 

In practice, the spectra of LRG contain several features that may be 
approximately modelled as breaks and/or line features of different wavelength widths. 
In the case of a single template, the optimal number of filters results from the 
contribution of all the features to the determination of the photometric redshift. At 
the same time, the spectra of real galaxies are, of course, not described by a single 
template and depend on a variety of parameters. The need to fit for these spectral 
variations in addition to the redshift favors an increased number of filters when the 
signal-to-noise that can be reached is high enough. This is seen in more detail in 
the following section; the analytic considerations discussed above will be helpful to 
interpret the numerical results for fits to simulated photometry of LRG. 

5.2. Single template galaxies 

We now present the numerical results of recovering photometric redshifts from the 
simulated photometry of LRG spectra, for the case of our fiducial survey as described 
in §4. The distribution of photometric redshift errors is generally not Gaussian because 
of the presence of outliers, or catastrophic failures, which naturally increase as the 
signal-to-noise drops. To avoid having to deal with identification of outliers, we report 
the photo-z error Cz in all of our figures, unless otherwise stated, as the interval such 
that 68% of the simulated galaxies have a smaller difference between their true redshift 
and fitted redshift. As long as the fraction of outliers is small, the 68 percentile value 
should approximately correspond to the rms redshift error after the outliers have been 
successfully removed. We shall not discuss here the methods for removing outliers 
in photometric redshifts; these would depend on how our assumed narrow-band filter 
system is complemented with wider filters on a broader range of wavelengths. In 
practice, any method for fitting photometric redshifts will fail to detect some outliers 
and will classify some good redshifts as outliers. Our assumption that outliers are 
removed perfectly is the most optimistic possible one. 

We start with the simplest case where the same template used to generate the 
galaxy photometry is then used to fit the photometric redshift. This corresponds to the 
ideal case where the galaxy spectral shape is precisely known, and only the redshift and 
luminosity need to be fitted. We use template 1 for this purpose in this section. The 
relative redshift error aiz = (7^/(1 + z) is shown in Figure [t] (top panel) as a function 
of redshift, for galaxies of luminosity and varying the number of filters. For this 
Figure, we show both the 68 percentile error {thick lines) and the 95 percentile error 
{thin lines). The thin and thick lines should be separated by a factor 2 for a Gaussian 
distribution of errors, and by an increasing factor as the fraction of outliers increases. 
The redshift errors are calculated at 450 values of the redshift, from z = 0.1toz = l 
in increments of Sz = 0.002. For each value of the redshift, the error is computed from 
a total of 100 photometry simulations. Therefore, each curve in Figure [t] (top panel) 
is based on a total of 45000 simulations of galaxy photometry. The results are shown 
as filled dots, after averaging in quadrature over redshift intervals of width Az = 0.1 
(lines between dots are plotted for guidance only) . The bottom panel of Figure [t] 
shows the redshift error (68 percentile) as a function of galaxy luminosity for redshifts 
0.65, 0.75 and 0.85 (with errors increasing with redshift), again after averaging over 
an interval Az = 0.1, and for Nf — 10, 30 and 50. 

We see from this Figure that for galaxies of luminosity L*, the redshift errors are 
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Figure 7. Top: Photo-z error as a function of redshift of an L* galaxy in our 
fiducial survey, varying the number of filters as indicated. The 68 or 95 percentile 
error is computed in each redshift bin of width Az = 0.002, and then averaged in 
quadrature over Az = 0.1. Bottom: Photo-z error (68 percentile, computed as in 
the top panel) as a function of galaxy luminosity, for redshift intervals centered at 
z = 0.65, 0.75 and 0.85 (errors increase with redshift), and for A^^ = 10, 30 and 50. 
Template 1 is used in all cases for generating and fitting the galaxy photometry. 
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practically independent of the number of filters in our fiducial survey. At z > 0.4, 
there is only a small improvement of cr^ when Nf is increased from 10 to 30, and 
even this small improvement is ofi^set by an increasing fraction of outliers as Nf is 
increased. This implies that it would not be worth to use more than ~ 10 filters to 
measure photometric redshifts for galaxies of luminosities near that are well fitted 
by a single template. 

The target accuracy of 0.003(1 + z) to measure radial BAO is achieved up to 
z ~ 0.7 for L = L^. The bottom panel shows that at 2; = 0.9 the same target can be 
achieved for L > 1.6L*. For galaxies of high luminosity, increasing Nf from 10 to 30 
becomes a stronger advantage, but gains for Nf > 30 are not substantial, especially 
because the outliers increase with Nf. The approximately constant redshift accuracy 
with the number of filters is roughly a consequence of the analytic arguments discussed 
in §5.1 for the case when the main spectral feature is an unresolved break. Note that 
there are also numerous features in the spectrum of template 1 (Figure 2) that can 
be modelled as absorption hues of width ~ A/30 (particularly at wavelengths just 
blueward of the Ha break), which are the principal reason for the improvement of 
redshift errors from Nf = 10 to Nf = 30 for high signal-to-noise. 

So far, the results have been shown averaged over redshift bins of Az — 0.1. The 
unbinned results are shown in Figure [8] for the cases Nf = 10 and Nf — 50. There 
are periodic oscillations in the redshift error, with a spacing Sz that corresponds to 
the separation between filters. These oscillations are due to different features in the 
spectrum (mostly the Ha break) crossing the center and the edge of filters as the 
redshift varies. The redshift error is best (smallest) when the Ha break is placed 
between two filters, and it is worst when it is placed at the center of one filter. These 
oscillations should naturally depend on the shape of the filter windows, as well as the 
form of the Ha break, and therefore the type of LRG galaxy. One should note that 
the scale of the oscillation corresponds to the BAO scale for Nf ~ 30. The presence 
of these artificial oscillations in the redshift error, which could vary their amplitude 
over the survey area owing to variations of the seeing or other observing conditions 
that may systematically affect the photometry, needs to be properly calibrated and 
corrected for in order to avoid introducing perturbations in the BAO peak of the 
galaxy correlation function. 

5.3. Using multiple templates 

In reality LRG are not a perfectly uniform population described by a single spectrum. 
We now see the impact on the redshift errors when variations in the galaxy spectra 
are allowed for. We use the five templates described in §3, which have been optimized 
to provide a good representation of observed galaxy spectra [20l [22] . 

Figure |9] shows the photo-z error for spectra generated with template 1 and 
L = L^. The thin lines are the same as those shown in Figure |7] (top panel) , fitted 
with the same template 1. The thick lines are the result obtained when any linear 
combination of the five templates is allowed in the fit. As expected, the errors increase 
when we allow for spectral variability: the maximum photo-z accuracy for an ideal 
galaxy having exactly the same spectrum as template 1 is obtained when we know 
that the galaxy has the spectrum of template 1. When we do not know this, and 
spectral variability needs to be allowed for, the photo-z error is larger. When fitting 
with our five templates, the target accuracy of az = 0.003(1 -I- z) is achieved only for 
z < 0.63 at i = for Nf = 30 (which is roughly the ideal number of filters). Results 
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Figure 8. Photo-z error as a function of redsliift witliout smoottiing for 10 filters 
(top) and 50 filters {bottom). The oscillations have a periodicity determined by 
the filter width. 




Figure 9. Impact of using multiple templates on redshift errors. Thin lines: The 
mock galaxy photometry is generated using template 1 and fitted with the 
same template. Thick lines: The same mock catalogs generated with template 
1 are fitted with all five templates. The 68 percentile errors are computed as in 
Figure [t] The number of filters is varied as indicated. 



for L > will be shown below in Figure 13 



The increase of the redshift error resulting from fitting with five templates is 
substantially worse for Nj — 10. When LRG are not assumed to be a homogeneous 
population, using a larger number of filters becomes more advantageous because of the 
need to distinguish the type of galaxy spectrum when fitting the redshift. However, 
as found in the previous section, increasing further the number of filters to Nf — 50 is 
generally not worthy, except at rather high signal-to- noise, where the number density 
of sources is small in our fiducial survey. 

The results in Figure |9] are a lower limit to the increase of redshift errors caused 
by galaxy variability, because real galaxies are not perfectly matched by the set of five 
templates we use. We discuss this further in §5.5. 

5.4- Variation of redshift error with templates and effects of galaxy evolution 



Figure 10 {top panel) shows the photo-z errors when the galaxy photometry is 
generated for each one of the five templates, and fitted allowing a linear combination 
of all five templates, for L = L^, and Nf = 30. There are substantial variations 
of the redshift error with the template spectrum, implying that the redshift error 
depends on the type of LRG galaxies selected to measure BAO. Over the redshift 
range 0.5 < z < 0.9, which is the most promising for our purpose, template gives 
the smallest errors and templates 2 and 4 the largest ones. The long-dashed line shows 
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Figure 10. Top panel: Mock catalogs of L, galaxies are generated with one 
template for each line (0-4) and fitted with the five templates. We also show the 
case where the mock catalog is generated by random linear combinations of the 
five templates and fitted with the five templates. Bottom panel: Photo-z accuracy 
for spectra of single starbursts of fixed age and metallicity and no dust, compared 
with template 0. The accuracy improves with increasing age and metallicity, and 
worsens with even trace amounts of star formation, but depends also on details 
of the modeling of giant and horizontal branch stars. 



Photo-z optimizaton 



19 



the error when random hnear combinations of the five templates are used to generate 
the galaxy photometry. As expected, the error obtained is an average of the errors 
from the templates to 4. 

To understand the reasons why the errors depend on the template used, we plot in 
the bottom panel of Figure 10 the photo-z error for a set of templates computed for a 
stellar population of fixed age (11, 5 and 2 Gyr) and metallicity {Zq, 0.2Zq, or 2.6^0); 
a template with an exponentially decaying star formation rate with r = 0.6 Gyr at 
age 5 Gyr is also included. These set of templates have been generated with the model 
of Bruzual & Chariot (2003), and are used only in this figure for illustrative purpose. 
Smaller errors are generally obtained for stellar populations that are older and of higher 
metallicity. The presence of even a small amount of young stars (the model shown has 
a star formation rate that has decreased by a factor e~^/°-^ = 2.4 x 10"''^ since star 
formation started) increases the errors. There are, however, other factors affecting 
the photo-z errors in addition to the age and metallicity of the stellar population. 
In particular, the modeling of the giant and horizontal branch stars affects the part 
of the spectrum near the Ha break, and has a complex effect on the redshift errors. 
We note that, for z > 0.5, template has smaller errors than all other templates. 
The variations between templates to 4 are due to differences in the distribution of 
stellar ages, amount of dust absorption, and different stellar evolution models used. 
Let us recall that templates to 4 are built using the empirically derived Panter et 
al. (2007) star formation histories. Template uses the SPEED models [H], which 
have been empirically found to provide a better fit to observations of LRG spectra 
than other stellar population models. Despite having the oldest stellar population, 
template is bluer than template 1 in the rest-frame wavelength range of 4000 to 
9000 A because it has a less developed giant branch and has no dust extinction. 
Template 1 aims at reproducing the spectrum of an SO/Sa galaxy, and has some dust 
absorption. Templates 1 to 4 are built using Chariot & Bruzual 2008 models, which 
have been empirically shown to provide a better fit to Sa-Sc galaxies. Their larger 
errors compared to template are due to different star formation histories (which vary 
among the templates in a complex way) and the different stellar evolution model. We 
note also that the redshift error depends on the band used to normalize the templates 
to the same luminosity; here we have normalized them to the same B-band luminosity, 
which is used to measure the luminosity function \}3 . 

The variation of the errors with the assumed stellar population in each template 
is roughly associated with the decreasing amplitude of the Ha break for younger and 
less metal-rich galaxies, and with the addition of trace amounts of star formation 
(see §5.1, and eq. [?]). For the same reason, changes in the modeling of giant and 
horizontal branch stars that affect the change of the Ha break have an important 
impact on redshift errors. Despite the complex dependence of the photo-z errors on 
several properties of the stellar population, the expected evolution with redshift is 
clear: stellar populations in LRG should be younger, star formation rates should 
increase, and the LRG population should become less homogeneous, leading to 
increased redshift errors. For similar reasons, redshift errors should decrease with 
galaxy luminosity. There is growing observational evidence that the most massive 
galaxies contain the oldest stellar populations up to z ~ 1 — 2 and they have less 
than 1% of their present stellar mass formed at z < 1 [551 1301 EJ [H] . In high-density 
regions (i.e., galaxy clusters), massive systems ceased their star formation by z '--^ 3 
[5T| . While galaxies with stellar masses > 5 x 10"'^^ M0 (roughly corresponding to 2L^,) 
ended their star-formation activity by z ~ 2, less massive objects were still actively 
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Figure 11. Photometric redshift versus input redshift. Our sample is generated 
by redshifting the spectra of SDSS DR6 LRG for L = L* and our fiducial survey 
(see text), and fitting them with our five templates. Dotted lines show the 
±0.05(1 + z) photometric error range. All points outside the dotted lines are 
considered outliers that could eventually be removed. 

forming stars at z ~ 1 [33j . As a result, the redshift errors found here using SDSS 
spectra of LRG, which are at z < 0.4 and L > 3L^,, are a lower limit to the redshift 
errors that may be found at higher redshifts and lower luminosities. 

5.5. Effects of variability in real LRG: SDSS Spectra 

In practical observations, the spectra of galaxies for which photometric redshifts are 
obtained are never exactly the same as the templates used for the fits. Real galaxy 
spectra are affected by the distribution in age, metallicity and dust obscuration, are 
not matched well enough by our imperfect models of stellar evolution, and may be 
substantially changed by very small amounts of star formation. For example, the 
presence of a weak Oil emission line (at a wavelength close to that of the Ha break) due 
to traces of star formation may change the recovered photo-z. In addition, superposed 
light from very faint galaxies along the line-of-sight may introduce small changes in the 
spectra. Therefore, the redshift errors we found so far by assuming that the adequate 
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Figure 12. Solid lines: Photo-z error for redshifted SDSS LRG at L = L« in 
our fiducial survey, fitted with our five templates. The thin line shows the 68% 
error in each redshift interval of width Az = 0.1 without removing any outliers, 
and the thick line shows the same error when the outliers outside the range of 
the dotted lines in Figure [TT] are removed. Dashed lines: Same results when we 
use a mock catalogue generated with template 1. The recovered photo-z errors 
for observed SDSS increase relative to the case of ideal galaxies generated exactly 
with one of the five templates used. 



templates are perfectly known may be severely underestimated. 

The best and most extensive available survey of LRG is the SDSS. In this Section 
we use observed spectra of LRG obtained from the publicly available SDSS-DR6; 
spectro . princeton . edu catalog. We randomly select 50, 000 galaxies flagged as LRG 
for which a high signal-to-noise spectrum is available. Unfortunately, most of the SDSS 
LRG are at redshifts z < 0.4. To simulate the measurement of a photometric survey 
at higher redshifts, we shift all the observed SDSS spectra by a constant factor of 
/ — 530/380 — 1.395 to longer wavelength. Hence, the true redshifts of the galaxies, 
z, are all shifted to a redshift for the simulated galaxy z' = (1 + z)/ — 1. With this 
shift factor, the minimum observed wavelength in the SDSS spectra, at A = 380 nm, is 
shifted just to the blue edge of our first filter at 530 nm, and the redshift range shifts 
to the interval 0.4 < z < 1.0, the most interesting range for measuring BAO with 
LRG photometric redshifts. We then follow the same procedure described in §4: the 
shifted spectra are convolved with the filter window functions, noise is added to the 
fluxes at each fllter, and the simulated photometry is fltted with our flve templates. 

We flrst examine a scatter plot of the fltted photometric redshift versus the true 
redshift of the SDSS galaxies (after the constant shifting mentioned above) for L — L^, 
and our flducial survey with Nf = 30, in Figure 11 For guidance, a dashed line for 
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z — Zphot has been added. The two dotted lines indicate 1 + Zphot = (1 ± 0.05)(1 + z), 
and will be our threshold for considering a photometric redshift as an outlier. The 
effect of the oscillations in the redshift error distribution with a period equal to the 
filter width, mentioned in §5.2, is also seen here. The dotted lines clearly separate 
a group of outliers in the error distribution, which arise because the Ha break is 
confused with other wider features in the spectrum when the noise is high. In the rest 
of this section, the results for the 68 percentile of the redshift error distribution will be 
shown for two different cases: considering the whole redshift error distribution, and 
removing the outliers defined as the points outside of the two dotted lines. The errors 
computed after removing the outliers should be considered as a highly optimistic case 
where the outliers are assumed to be perfectly identified with the help of additional 
broad-band filters in the survey. 

The solid lines in Figure [12] are the 68 percentile redshift error as a function of 



redshift, for L — and our fiducial survey with Nf = 30, as in Figure 11 The thin 
line is the value for the whole distribution, and the thick line is obtained after the 
outliers are removed. The effect of the outliers is substantial, increasing the redshift 
error by ~ 40%. For comparison, we show as dashed lines the case where galaxies 
generated with template 1 are fitted allowing for a linear combination of all our five 
templates, also for the cases of removing the outliers {thick line) or not removing them 
{thin line). Note that the thin dashed line is not exactly the same one as the thick 
solid line in Figure [9] the reason is that here, the 68 percentile error is found over each 
redshift interval of width Az = 0.1, in the same way as for the SDSS galaxies, whereas 
in previous plots the 68 percentile error is found first in bins of width Az — 0.002 
and then averaged in quadrature over the larger width Az = 0.1. The errors are 
slightly reduced for this reason in Figure [12] and in reality one would have to deal 
in some way with the large oscillations with redshift shown in Figure |8] which are 
the main reason for the change in the error with the redshift bin width used for 
extracting the 68 percentile value. Even after the outliers are removed, the redshift 
errors obtained when fitting the realistic LRG population are substantially increased 
relative to fitting an ideal population matching template 1 exactly. An error below 
0.003(1 + z) is obtained only at 2; < 0.55 for L — L^,. The rms error is still larger. 



because even after removing the outliers outside the dotted lines of Figure 11 the 



error distribution is still substantially non-Gaussian. Figure [12] also shows that the 
number of outliers is much smaller for galaxies generated with template 1. Real galaxy 
variability introduces complications in fitting photometric redshifts that increase the 
fraction of outliers. 

Overall, the 68 percentile redshift error for galaxies increases by about 20% 
for SDSS galaxies compared to template 1 galaxies after the outliers are removed, 
and by ~ 40% if the outliers are not removed. Comparing also to the results when 
template 1 galaxies are fitted with the same template 1 in Figure [9] we find that the 
redshift errors for the SDSS galaxies are 50% to 60% larger than the errors found 
when assuming that there is no galaxy variability at all, after removing the outliers, 
and a factor of 2 larger when outliers are not removed. 

The redshift errors after removing outliers are also shown in Figure |13| as a 
function of galaxy luminosity, for several redshifts. Results are presented for the 
SDSS galaxies and for the cases of fitting template 1 galaxies with either template 1 
alone or all our five templates. The error in the SDSS galaxies decreases much more 
slowly with luminosity than for mock galaxies generated with the templates. This is 
simply because the variability of SDSS LRG spectra is not well represented by our 




Figure 13. Photo-z error as a function of luminosity. Solid lines: Redshifted 
SDSS DR6 galaxies are fitted with our five templates. Dashed lines: Same results 
when galaxies are generated with template 1, and fitted with five templates. 
Dotted lines: Same results for galaxies generated with template 1 and fitted 
with the same template. The redshift error is the 68 percentile of the averaged 
distribution over each redshift interval of width Az = 0.1, centered at z = 0.65, 
0.75 and 0.85 (errors increase with redshift). 



templates. For L > 2L^, the errors for the SDSS galaxies are nearly twice as large as 
the errors obtained from galaxies that match template 1 exactly, while the errors when 
fitting with one or all five templates become nearly equal. The fraction of outliers is 
also found to be much larger for the SDSS galaxies. The galaxy variability in the 
SDSS LRG is highly complex and cannot be easily accounted for with a small number 
of templates, even when these templates have been calibrated with observed galaxies. 

We have checked the dependence of the redshift errors on the signal-to-noise 
reached in the photometry. For Nf — 30, we find that the photo-z error versus the 
signal-to-noise in the reddest filter {S/Niast) is quite robust to the various assumptions 
made about the survey design. The target accuracy of ai^ = 0.003 is reached for SDSS 
galaxies at S/Niast — 20, and at S/Niast — 12 in the one template anlaysis. The photo- 
z error and the fraction of outliers increase rapidly for lower signal-to-noise. At higher 
signal-to-noise, the errors continue to decrease rapidly for galaxies generated with the 
templates, but they flatten out to an "error-floor" of about cr/^ = 0.002 for the SDSS 
galaxies, the "knee" being at around S/Niast ~ 30, indicating an intrinsic galaxy 
variability that is not adequately modelled by a small number of templates. 

A possible problem when using the SDSS spectra of LRG to model photometry 
might arise from spectroscopic calibration errors, which could result in artificial 
variations of the simulated fluxes that are not due to real galaxy variability but 



Photo-z optimizaton 



24 



to calibration errors. In fact, [16] argued that SDSS spectra cannot be used to 
model variability for this reason; we reach a different conclusion. Even though 
the spectroscopic calibration errors are found to be typically at the level of ^ 4% 
on wavelength scales of 100 nm when comparing photometry and spectroscopy of 
stars (see |http://www.sdss.org /dr6/pro ducts/spectra/spectrophotometry.html| ), the 
precision of the photometric redshifts depends on measuring spectral features at the 
smallest wavelength scale allowed by the narrow-band filters (for LRG, this is mostly 
the Ha break). For a filter width of lOnm, the derived photometric redshifts might 
be affected by spectroscopic calibration errors on a similar wavelength scale. At this 
scale, the calibration errors are much smaller: they are certainly less than 2% from 
the analysis of quasar spectra in regions with no Lya forest absorption (P. McDonald, 
private communication; see also Figure 19 in [32], where the contribution to the flux 
variance due to calibration errors from different scales in quasar spectra is shown 
as a ratio to the Lya forest power), and they should be less than 1% from the 
most recent calibration made with low-metallicity F-subdwarfs (D. Schlegel, private 
communication). We therefore conclude that the variations in the SDSS spectra of 
LRG galaxies which result in increased photometric redshift errors are mostly real and 
not due to spectroscopic calibration errors. 

We emphasize that, owing to galaxy evolution, the variability present in the 
SDSS spectra is a lower limit for the variability that should be expected for LRG at 
higher redshift, as discussed in §5.4. The errors might possibly be reduced by further 
optimization of the templates as a function of galaxy luminosity and environment if 
galaxy variability at high redshift were better understood, but how substantial an 
improvement may be achievable in this way remains to be proven. 



6. Discussion 



This paper has presented an analysis of the expected photometric redshift errors in 
LRG that can be achieved in an imaging survey with a system of a large number of 
narrow-band filters. We have found that an optimal choice for the number of filters 
is Nf ~ 30 over the wavelength range 530 to 830 nm, allowing redshift measurements 
over the range 0.5 < z < 0.9. A central focus of our investigation is that the intrinsic 
variability of LRG implies a substantial increase of the redshift errors compared to 
a calculation where all the LRG are assumed to match exactly the same template 
used to fit their photometry. We find that it is crucial to properly include the effects 
of spectral variability when evaluating the prospects for any survey to measure BAO 
(particularly in the radial direction) in the power spectrum of LRG. 

The implications of our results can be summarized by considering, for any 
particular survey, the lowest luminosity that a galaxy needs to have to yield a redshift 
accuracy better than the target level for measuring BAO in the radial direction, 
aiz < 0.003. The power of a survey to reach the required photometric accuracy 
depends mainly on the product Et of the etendue times the total observing time. 
This lowest luminosity as a function of Et is shown in the upper panel of Figure 14 at 
five different redshift intervals, when the total survey area is fixed to = 8000 square 
degrees and Nf = 30, for the SDSS galaxies after removing outliers with redshift errors 
Az > 0.05(1 -I- z). Note that the minimum luminosity shown at each value of Et is the 
one that would yield a 68 percentile error equal to aiz = 0.003; the rms error would be 
larger because even after removing the most extreme outliers, the error distribution 
is still substantially non-Gaussian. The luminosity function [19] implies the number 





Figure 14. Top panel: Luminosity of an LRG with photometric redshift error 
ai^ = 0.003 after removing outliers, as a function of the etendue X observing 
time of the survey, for the five redshift intervals indicated. Errors are computed 
using the spectra of SDSS LRG. Bottom panel: Corresponding number density of 
galaxies with a photometric redshift error better than 0.003(1 + z), as a function 
of Et, from the luminosity function |19| after taking out the fraction of outliers. 
Horizontal dashed lines indicate the required galaxy densities for reaching a shot 
noise with nP = 1 at two wavenumbers (differing by a factor of two in the total 
number of modes they comprise), where P is the linc-of-sight power spectrum 
assuming ai^ = 0.003(1 + z) and shown in Figure ^ Our fiducial survey has 
Et = 1.5 X 10^ m2 deg2 hr. 
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density of galaxies above this minimum luminosity that is shown in the bottom panel, 
after subtracting outliers. As shown in Figure [T] the photometric redshift errors 
suppress the power spectrum near the line-of-sight, which must be compensated by an 
increased number density of galaxies to keep shot noise sub-dominant on a given scale. 
The bottom panel of Figure 14 addresses this issue. For redshift errors aiz = 0.003, 
the two horizontal dashed lines show the galaxy number density required to reduce 
the shot noise to the level nP = 1 at wavenumbers k = 0.13 h Mpc~^ and k = 0.164 
h Mpc~^ near the line-of-sight. Modes up to fc = 0.164 h Mpc~^ must be sampled 
to fully measure the second BAO. Figure 14 shows that Et > 3 x 10^ m^ deg^ hr 
is required to reach nP = 1 up to z = 0.85 for k = 0.164 h Mpc^^ near the line-of- 
sight. For the characteristic total observing time of 5000 hours assumed in our fiducial 
survey, this corresponds to i!^ > 60 m^ deg^. 

This minimum value for Et depends on several technical details of the telescope- 
camera system and the observing conditions for the survey. The assumptions made 
here for our fiducial survey were specified in §4: a seeing of 0.8", read-out noise 
of 7 X 3 electrons per pixel (with a pixel size of 0.4"), the efficiency for the latest 
available CCD's, optical losses of two mirror reflections only, and the sky brightness 
of the Paranal observatory at an airmass of 1.2. These assumptions can be considered 
as highly optimistic, except perhaps the assumption on the read-out noise which might 
be improved with the best available technology. Changes in these assumptions imply 
corresponding changes in the required value of Et necessary to reach a fixed signal- 
to-noise in the galaxy photometry. We can quantify this in the following way: 

• A seeing degradation from 0.8" to 1.2" implies a 20% decrease in the signal- 
to-noise for a fixed aperture, and therefore requires a 40% increase in Et 
{din Et / dhi seeing" ~ 1) 

• An increase by half a magnitude of the sky brightness (comparable for example 
to the sky variation over a sun cycle or between ecliptic plane and ecliptic pole) 
can be compensated by a increase in Et of a factor 1.6 {din Et/dBsky — 0.9) 

• A reduction in readout noise of a factor of three can be compensated by a ~ 20% 
reduction in Et { din Et/ din N readout — 0.3; we have verified this by repeating all 
the above calculations for a lower read-out noise). An even lower read-out noise, 
however, would not yield the same fractional improvement because at that point 
the read-out noise becomes subdominant. 

• Any factors that decrease the overall throughput of the system (i.e., lower 
efficiency of the CCD's, additional optical losses in the telescope-camera system, 
reduced fllter transmission, etc.) will of course increase the required Et in 
proportion to the inverse of the throughput. An important factor that may cause 
such a reduction in the throughput would likely be present in a system of narrow 
filters based on interferometry, where the filter window has a variation across the 
field-of-view caused by the varying incidence angle of the light. The averaging 
of the filter window shape as a source moves across the field in drift-scanning 
mode would widen the filter window without increasing the number of detected 
photons. 

These considerations mean that, in practice, the required etendue would likely be 
substantially larger. 
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7. Conclusions 

The principal objective of this paper is to examine the potential for narrow-band 
photometric surveys to measure the radial BAO signal in the LRG power spectrum. 
We have revisited this issue after [T5J [TTJ [T3], addressing the optimization of the 
number of filters and filter width, discussing the dependence on the spectra of 
the target population of LRG, and considering the signal-to-noise that could be 
realistically achieved in a survey. We have found it to be particularly important to 
take into account the spectral variability of a realistic galaxy population for properly 
evaluating the photometric redshift accuracy that can be achieved. Our conclusions 
can be summarized as follows: 

a) In agreement with |12lllllfn] the photometric approach, even if it can achieve 
the target photo-z error of 0.3%, will be advantegeous only if it can cover a much 
larger fraction of the sky, with a higher galaxy density and reach higher redshifts 
than the spectroscopic surveys currently under way. Not only one needs higher source 
density in a photometric survey to measure the oscillations up to fc ~ 0.2, but also 
the distribution of redshift errors needs to be known accurately to correct for their 
effect on the power spectrum shape and be able to measure the BAO shape. This is 
equivalent to the requirement of obtaining photo-z with better than 0.3% accuracy 
to fainter magnitudes than spectroscopy. We have concentrated on the LRG galaxies 
because they are a fairly homogeneous population, the form of their spectra makes 
them particularly suitable for good photo-z determinations and because they represent 
the target of choice for z < 2 BAO surveys. 

b) If LRG were a perfectly homogeneous population (with spectra closely matched 
by a single spectral template), the photometric redshift accuracy would not improve 
very much beyond a total number of filters Nf ~ 10, except for galaxies at high 
signal-to-noise which would have a low number density in a survey with Et ~ 
1.5 X lO^m^deg^ hr. In reality, galaxy spectra are variable and the ideal number 
of filters is Nf ~ 30 over a wavelength range 530 nm < A < 830 nm, focusing on 
the redshift range 0.5 < z < 0.9. Medium-to-narrow band filters photometric redshift 
errors show periodic oscillations with a spacing Sz that corresponds to the filter widths. 
For a filter width of ^ 100 A, the scale of the oscillation corresponds to the BAO scale. 
The presence of these oscillations may introduce significant biases in the measurement 
of the BAO scale (see Section 5.2). 

c) The variability of realistic galaxy spectra demand high resolution (or more 
filters in a photometric survey) and high signal-to-noise. We have quantified, with 
realistic simulations of a photometric survey, the minimum galaxy luminosity and 
the corresponding number density for which 0.3% photo-z error can be obtained as a 



function of the survey Etenduex exposure time Et (see Figure 14 1. We estimate the 
minimum value of Et necessary to measure the full second baryonic acoustic oscillation 
in the power spectrum, at fc < 0.164 h/Mpc near the line-of- sight, with shot noise 
nP > 1 and up to z = 0.85, at Et ~ 3 x 10^ m^ deg^ hr. The actual required value of 
Et would likely be larger when taking into account that our assumptions for the seeing 
and the overall throughput of the telescope-camera system are optimistic, and that 
the presence of outliers, optimistically clipped and ignored, will degrade the estimated 
performance. Moreover, the true level of galaxy variability should likely be higher than 
in the SDSS LRG sample we have used (which is for z < 0.4 and L > 3L*), because 
the variability of LRG spectra is expected to increase with redshift and decrease with 
luminosity. 
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Current spectroscopic surveys (e.g., SDSSIII-BOSS) are obtaining spectra of the 
massive (> 2L*) population of LRG at z < 0.75 to measure radial and tangential 
BAO scale. Here we have investigated the conditions necessary to obtain the required 
photo-z accuracy of 0.3% from a similar or fainter population using photometry, which 
in principle is less expensive in terms of observing. We have concluded that an etendue 
of 30 m^ deg^ is required for tracing the same population as BOSS {z < 0.75, L > 2L*) 
with radial BAO signal degraded only by 60% due to photo-z errors. But to improve 
on this i.e., to reach galaxies at slightly higher redshifts, the minimum etendue 
required is 60 m^ deg^ for a 5- year survey (corresponding to a 4m telescope with the 
largest fields of view that have been made), and probably needs to be larger given 
our optimistic assumptions we have made for the seeing, throughput, observing time, 
removal of outliers and systematic errors. 

In this paper we have restricted our attention to the objective of measuring radial 
BAO with LRG. However, a large-area imaging optical survey with a large number of 
narrower bands than existing surveys would have many other applications, and might 
lead to a large number of interesting astronomical discoveries. We have shown in this 
paper that, for realistically achievable values of Et at the present time, the density of 
LRG sources that could be measured with a narrow-band photometric survey with the 
target redshift accuracy for radial BAO will not be large enough to make it competitive 
with a spectroscopic survey. Therefore, a photometric survey with a large number of 
optical narrow-bands needs to find its scientific justification in other astronomical 
applications. 
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